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LONG-TERM  GOALS 

The  overall  goal  of  this  project  is  to  determine  the  feasibility  of  in  vivo,  non- invasive  measurement  of 
the  complex  elastic  moduli  of  cetacean  head  tissue.  If  this  objective  is  met,  measurement  systems 
could  be  developed  capable  of  measuring  the  complex  elastic  moduli  of  the  head  tissue  of  live, 
stranded  cetaceans. 

OBJECTIVES 

The  technical  objective  is  to  remotely  generate  and  detect  mid- frequency  (I  to  lOkHz)  elastic  waves 
within  the  body  of  a  living  cetacean,  using  ultrasound  and  to  use  the  measured  propagation  parameters 
of  these  waves  to  obtain  the  complex  elastic  moduli  by  inversion.  A  further  technical  objective  is  to 
extract  tissue  moduli  in  this  manner  intracranially.  This  objective  carries  considerably  more  technical 
risk  since  both  the  wave-generating  ultrasound  and  the  probe  ultrasound  will  be  attenuated,  distorted 
and  scattered  by  the  passage  through  the  skull. 

APPROACH 

The  approach  is  to  measure  the  complex  shear  and  bulk  modulus,  from  which  all  other  moduli  can  be 
calculated.  The  shear  modulus  will  be  determined  by  measuring  the  speed  and  attenuation  of  shear 
waves  generated  within  the  tissue  using  focused  ultrasound  as  a  remote  localized  force  generator. 

This  general  approach  to  determining  the  complex  modulus  is  an  application  of  a  new  medical  imaging 
technology  called  elastography.  The  methods  described  by  Greenleaf  (Chen  et  al,  2002)  and 
Sarvazyan  (Sarvazyan  et  al,  1998)  provide  the  basis  for  shear  wave  generation.  Displacements  are 
generated  remotely  in  a  tissue  volume  using  one  or  more  focused  ultrasound  beams.  In  the  single¬ 
beam  case,  the  ultrasonic  carrier  signal  is  modulated  at  a  low  frequency  Af.  In  the  dual-beam  confocal 
configuration,  the  two  ultrasonic  signal  carrier  frequencies  are  offset  from  each  other  by  2Af.  In  both 
cases,  a  radiation  force  at  the  focal  point  primarily  generates  shear  waves  at  the  frequency  2Af. 
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The  particle  displacements  resulting  from  the  passage  of  the  propagating  shear  wave  can  be  detected 
remotely  using  a  modified  version  of  an  ultrasonic  Doppler  vibration  measurement  system  called 
NIVAMS  developed  at  Georgia  Tech  (see  Cox  and  Rogers,  1987,  Martin  et  al,  2002).  The  system  will 
be  modified  to  make  the  transmitting  and  receiving  beams  coaxial  in  order  to  permit  operation  in  tight 
quarters.  Algorithms  will  be  developed  to  enable  the  magnitude  and  phase  of  vibration  to  be 
determined,  as  well  as  the  range  at  which  the  vibration  is  being  measured.  By  measuring  the  amplitude 
and  arrival  time  of  the  shear  wave  at  two  different  points  the  propagation  speed  and  attenuation  can  be 
determined. 

Elastic  waves  will  be  both  remotely  and  directly  generated  in  tissue  phantoms  and  measured  both 
remotely  and  directly  to  validate  the  measurement  technique.  The  elastic  properties  of  tissue  phantoms 
will  be  obtained  from  remotely  generated  and  measured  data  and  compared  with  directly  measured  and 
tabulated  material  values.  The  noninvasive  technique  will  be  repeated  for  tissue  phantoms  enclosed  in 
a  simulated  or  hydrated  real  cetacean  skull,  and  with  harvested  tissue  samples.  Ultimately,  feasibility 
will  be  assessed  through  measurements  on  a  Navy  bottlenose  dolphin.  The  animal  tests  will  be  done 
initially  for  accessible  tissue  and  subsequently  for  intracranial  tissue.  Ultrasound  intensity  and 
exposure  time  will  be  consistent  with  limits  that  have  been  established  as  safe  for  humans  and 
frequencies  will  be  kept  high  enough  to  be  far  above  the  highest  frequency  that  is  audible  to  the 
animal.  The  measurements  should  be  harmless  and  painless. 

WORK  COMPLETED 

1.  Ultrasonic  Vibrometer  Development  A  prototype  range-resolving  ultrasonic  vibrometer  was 
developed  and  its  performance  was  demonstrated  through  a  series  of  tests  with  a  tissue  phantom  which 
mimics  the  ultrasonic  properties  of  mammalian  brain.  The  system  concept  was  validated  through  the 
observation  of  shear  wave  displacements  and  calculation  of  their  propagation  speed  from 
measurements  of  delay  as  a  function  of  position. 

2,  System  Characterization  The  behavior  of  the  combined  ultrasonic  force  generation  and  vibration 
measurement  systems  was  examined  in  the  laboratory  using  a  suspended  point  scatterer.  This  was 
particularly  useful  for  characterizing  the  force  generation  system  in  terms  of  drive  level  and  effective 
beamwidth. 

RESULTS 

1.  Ultrasonic  Vibrometer  Development 

In  previous  work,  two  different  techniques  have  been  used  for  ultrasonic  vibrometry.  The  first  of  these 
involved  a  pure  tone  ultrasonic  carrier  signal  from  which  a  phase  modulation  produced  by  the 
displacement  of  a  reflecting  surface  was  extracted  to  compute  the  time  history  of  the  surface 
displacement.  This  technique  permitted  the  transduction  of  very  small  displacements  (amplitudes  on 
the  order  of  a  millionth  of  the  ultrasonic  wavelength)  but  did  not  allow  for  differentiation  between 
individual  scatterers  contributing  to  the  received  signal.  Thus  it  would  not  have  permitted  the 
interrogation  of  a  weak  scatterer  such  as  cetacean  brain  tissue  directly  behind  a  strong  scatterer  such  as 
bone.  The  second  technique  involved  a  pulsed  ultrasonic  carrier  signal  from  which  the  motion  of  a 
reflecting  surface  was  extracted  by  repeated  correlation  of  the  echoed  pulses  with  those  that  were 
transmitted.  This  technique  permitted  different  scatterers  to  be  simultaneously  interrogated,  but 
required  that  the  motion  of  the  scatterers  be  relatively  large  (amplitudes  great  than  a  tenth  of  a 
wavelength).  Shear  motions  of  this  size  cannot  be  practically  or  safely  generated  in  vivo.  The 
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technique  that  is  currently  being  developed  is  intended  to  combine  the  measurement  resolution  of  pure- 
tone  vibrometry  with  the  spatial  resolution  of  pulse-echo  vibrometry. 


The  new  technique  requires  that  the  ultrasonic  carrier  be  a  continuous  signal  made  up  of  a  large 
number  of  discrete  tones  that  have  a  high  spectral  purity  (less  than  -100  dB  spectral  side-lobes  at 
greater  than  2  Hz  separation  from  each  tone).  This  allows  for  the  received  signal  to  be  decomposed 
into  two  components  with  different  time  dependences  by  using  a  novel  spectral  interpolation  and  pulse 
compression  technique.  One  of  these  components  (the  carrier)  is  dependent  only  on  the  high-frequency 
time  scale  of  the  ultrasonic  signal,  and  the  other  (the  modulation)  is  dependent  only  on  the  low- 
frequency  time  scale  of  the  acoustic  path  length  modulation.  The  two  components  can  be  separated  in 
the  frequency  domain  and  then  be  recombined  to  give  a  displacement  history  of  all  the  scatterers 
contributing  to  the  received  signal  as  a  function  of  the  mean  acoustic  path  length  to  each.  With 
appropriate  hardware  and  signal  post-processing,  this  technique  yields  displacement  resolution 
comparable  to  pure-tone  vibrometry  at  the  center  frequency  of  the  carrier  signal  and  spatial  resolution 
superior  to  pulsed  vibrometry  exploiting  a  comparable  overall  bandwidth.  This  improvement  in  spatial 
resolution  is  achieved  by  combined  the  results  of  two  pulse  compressions,  which  make  use  of 
compression  waveforms  that  are  in  quadrature  with  each  other.  The  functional  representation  of  this 
combination  is  as  follows: 
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Where  5^{t\^  is  the  displacement  of  the  scatterer,  the  /  and  Q  values  represent  the  result  of  the 

two  pulse  compressions,  the  car  and  mod  subscripts  indicate  the  separated  carrier  and  modulation 
terms  respectively,  and  Ch  is  the  propagation  speed  of  the  ultrasonic  signal.  The  relationship  is  valid  for 
th  corresponding  to  the  two-way  propagation  time  to  the  scatterer.  Implementation  of  equation  1 
gives  a  complete  time  history  of  the  displacement  in  pseudo  real  time  (processing  time  is  currently 
about  ten  times  data  acquisition  time  on  a  desk  top  PC,  but  this  is  limited  primarily  by  memory  access 
speed). 

As  implemented,  the  prototype  system  exploited  a  carrier  signal  with  101  discrete  tones  separated  by 
6.25  kHz  and  centered  on  2.5  MHz.  This  sampled  low-frequency  vibration  in  the  range  of  2  Hz  to  3.12 
kHz  at  960  points  within  a  brain-mimicking  tissue  phantom  (Blue  Phantom  Inc.).  The  measurement 
points  spanned  from  the  face  of  concentric  transmitting  and  receiving  transducers  to  a  depth  of  1 1.5  cm 
inside  the  phantom.  The  measurement  depth,  measurement  resolution,  and  the  ability  to  separately 
depict  motion  at  each  individual  point  depend  on  a  variety  of  factors  related  to  properties  of  both  the 
tissue  and  the  transmit  signal.  Attenuation  must  be  both  sufficiently  small  to  permit  measurements  to 
the  full  depth  possible  with  the  tonal  spacing  of  the  transmitted  signal  and  sufficiently  large  to  preclude 
spatial  ambiguities  caused  by  wrap  around  phenomena  in  the  integral  transforms  used  in  the  pulse 
compressions.  This  in  turn  dictates  the  center  frequency  of  the  ultrasonic  carrier  signal  as  attenuation  is 
highly  frequency  dependent.  The  ability  to  separate  points  in  space  depends  on  homogeneity  (speckle 
scatter)  of  the  tissue  and  on  the  bandwidth  of  the  carrier.  In  general,  an  increased  bandwidth  improves 
the  spatial  resolution  by  permitting  compression  pulses  with  shorter  duration.  In  practice,  this 
bandwidth  is  limited  by  hardware.  It  is  difficult  to  distinguish  the  motion  of  a  region  of  weak  scattering 
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in  close  proximity  to  a  region  of  strong  seattering  in  mueh  the  same  way  that  these  regions  are 
shadowed  in  conventional  pulse-eeho  ultrasonic  imaging. 

Initial  testing  of  the  ultrasonie  vibrometer  involved  direet  meehanieal  exeitation  of  shear  waves  in  a 
tissue  phantom  suspended  in  a  water  filled  tank,  as  shown  in  Figure  1 .  This  was  done  in  order  to 
simplify  the  experimental  system  and  faeilitate  refinement  of  the  hardware  and  processing  algorithms. 
In  this  eonfiguration,  the  phantom  was  driven  on  one  side  by  a  19mm  polycarbonate  sphere  that  was 
compelled  to  an  eleetrodynamic  shaker  by  means  of  a  narrow  stinger.  The  motion  of  the  sphere  was 
along  an  axis  at  a  shallow  angle  to  the  surfaee  of  eontaet.  This  exeited  an  outwardly  propagating 
spherieal  shear  wave  in  the  phantom  with  an  amplitude  that  was  relatively  constant  over  a  large 
volume  and  eould  be  eontrolled  to  mateh  the  amplitudes  expeeted  from  radiation  force  excitation.  The 
vibrometer  was  implemented  with  a  dual-ring  eonfoeal  ultrasonie  transdueer,  with  one  ring  radiating 
the  earrier  signal  and  the  other  ring  reeeiving  the  backseattered  pressure. 


Figure  1.  Vibrometer  test  setup.  Shear  waves  were  generated  by  an  external  vibration  source  (S). 

The  vibrometer  transducer  (T)  was  attached  to  a  positioner  (P)  to  facilitate  scanning.  The  phantom 

(B)  and  transducer  were  submerged  in  a  water  tank. 

Figure  2  shows  the  displacements  from  one  of  these  measurements  at  all  960  points  plotted  in  groups 
of  120  on  a  ±4pm  seale.  Here  the  shear  exeitation  signal  was  a  swept-frequeney  chirp  and  the  data  was 
analyzed  in  a  band  from  2  Hz  to  200  Hz.  Close  to  the  ultrasonie  transdueers  (upper  left  of  the  figure) 
the  measurement  is  overwhelmed  by  noise.  This  region  is  in  the  water  outside  the  tissue  phantom.  The 
high  noise  floor  here  is  a  eonsequence  of  weak  scattering  from  the  water.  Mathematically,  this 
corresponds  to  a  null  in  the  denominator  of  equation  1 .  Over  the  remainder  of  the  beam  the  measured 
time  waveform  elearly  eorresponds  to  the  exeitation  signal.  The  interferenee  midway  through  the 
waveform  is  attributable  to  the  shear  wave  refieetion  from  the  bottom  surfaee  of  the  phantom.  At  long 
ranges  the  signal-to-noise  ratio  diminishes  due  primarily  to  a  rising  noise  floor  rather  than  a  falling 
signal  level.  This  is  beeause  the  attenuation  in  the  phantom  reduces  the  energy  returned  in  the 
ultrasonic  carrier  at  longer  ranges. 
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Figure  2:  Measured  displacement  time  histories  at  960  locations  along  an  ultrasonic  interrogation 
beam  inside  a  tissue  phantom.  Signals  shown  in  the  upper  left  plot  are  in  the  water  path  between 
the  transducers  and  the  phantom.  Signals  shown  in  the  lower  right  plot  are  furthest  into  the 
phantom  and  indicate  increased  noise  due  to  attenuation  of  the  ultrasonic  carrier. 


In  Figure  3  the  delay  times  resulting  from  a  eross  eorrelation  of  eaeh  measured  displaeement  time 
waveform  with  the  shaker  drive  signal  are  plotted  for  all  of  the  loeations  along  the  beam.  Here,  the 
ultrasonie  transducers  were  located  three  different  distances  away  from  the  spherical  shear-wave 
source.  Also  depicted  on  this  figure  is  the  expected  delay  over  these  locations  that  was  computed  from 
a  separate  resonant-bar  measurement  of  the  bar  wave  speed  characterizing  the  tissue  phantom  material. 
The  measured  data  from  the  ultrasonic  vibrometer  can  be  seen  the  match  this  prediction  quite  well.  The 
observed  curvature  of  the  delay  measured  along  each  beam  individually  also  corresponds  well  to  the 
distances  from  the  source  to  each  measured  point.  This  is  because  the  wavefronts  were  spherical  and 
the  beam  linearly  transected  these. 

In  summary,  the  technique,  hardware,  and  processing  algorithms  for  a  high  precision  depth- 
discriminating  ultrasonic  vibrometer  have  been  developed  and  tested.  The  prototype  system  is  capable 
of  distinguishing  the  motion  of  regions  a  few  millimeters  apart  and  characterizing  that  motion  to  a 
precision  of  approximately  1  nm/VHz  with  a  2.5  MHz  center-frequency  carrier  signal.  The  practical 
limitations  of  this  system  are  well  understood  and  it  appears  to  be  suitable  to  the  in  vivo  interrogation 
of  a  cetacean  brain. 
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Figure  3:  Delay  times  as  a  function  of  depth  in  the  ultrasonic  interrogation  beam  for  three 
transducer  locations.  The  measured  propagation  delay  is  in  good  agreement  with  the  expected 
value.  Measurement  locations  less  than  approximately  1.0  inch  from  the  transducer  are  in  the 

water  path  between  the  transducer  and  phantom. 

2.  System  Characterization 

The  combined  radiation  force  generation  and  vibrometer  systems  were  characterized  using  the 
laboratory  setup  shown  in  Figure  4.  The  inner  transducer  was  used  as  the  vibrometer,  and  the  outer 
transducer  was  used  as  the  radiation  force  generator.  A  point  scattering  target  was  provided  by 
suspending  a  0.020”  diameter  stainless  steel  ball  on  a  synthetic  filiment.  Each  transducer  was  fixed  to 
a  positioner  to  allow  alignment  and  subsequent  field  scanning  by  translation  of  one  transducer  relative 
to  the  other.  This  is  the  general  approach  for  the  in  vivo  feasibility  assessment  prototype  system'. 

The  test  setup  was  used  to  check  a  variety  of  parameters,  including  force  field  shape  as  a  function  of 
drive  configuration  (amplitude  modulated  uniform  drive  vs  dual  ring  offset  carrier),  “linearity”  of 
induced  displacement  relative  to  the  square  of  drive  voltage,  minimum  detectible  displacement,  and 
measurement  repeatability.  Figure  5  shows  the  measured  displacement  of  target  as  the  forcing 
transducer,  operated  in  the  amplitude  modulation  configuration,  is  translated  laterally  away  from  the 
center  of  focus.  The  displacement  (“direct”)  is  normalized  to  the  peak  value  and  plotted  along  with  the 
expected  beam  pattern  (“calibration”)  based  on  a  separate  measurement  of  the  radiated  power  from  the 
transducer.  The  main  lobe  half-width  of  approximately  1mm  is  as  expected.  The  observed 
displacement  sidelobes  are  approximately  25%  higher  than  predicted,  although  this  is  not  expected  to 
materially  impact  system  performance.  The  data  collected  in  the  characterization  runs  showed  the 
ability  to  remotely  vibrate  targets  and  measure  the  resulting  displacements  with  an  all-ultrasonic 
system,  providing  a  basis  for  combined  system  testing  with  tissue  phantoms  now  underway. 


'  For  testing  on  tissue  samples  and  marine  mammals,  a  compact,  lightweight  positioner  has  been 
designed  and  built  so  that  the  device  could  be  hand-held. 
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(a)  (b) 

Figure  4.  System  characterization  test  setup,  (a)  nested  transducer  configuration,  (b)  cross-section 

illustration  with  transducers  focused  on  a  point  target 


Displacement  Beam  Pattern 


Figure  5.  Radiation  force-induced  displacement  beam  patterns.  “Direct”  refers  to  the  result 
obtained  with  the  ultrasonic  force  generation  and  reception  systems.  “Calibration  ”  refers  to  a 
prediction  based  on  a  separate  measurement  of  the  forcing  transducer’s  radiated  pressure  field.. 

IMPACT/APPLICATIONS  (Potential  future  impact  for  Science  and/or  Systems  Applications) 

There  is  considerable  interest  in  the  development  of  structural  acoustic  models  for  the  cetacean  head 
for  two  main  reasons:  1)  to  better  understand  biomechanics  of  sound  reception  and  production  in 
cetaceans,  and  2)  to  understand  and  hopefully  mitigate  any  harmful  effects  of  man-made  sound  on  their 
health  and  behavior.  The  development  and  validity  of  these  models  is  severely  limited  by  an  almost 
complete  lack  of  knowledge  of  the  mechanical  properties  of  the  constituent  tissue.  It  is  known  that 
these  properties  change  upon  the  death  of  the  animal.  There  is  thus  considerable  interest  in  being  able 
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to  measure  these  properties  in  vivo.  The  teehniques  and  instrumentation  investigated  here  should  also 
have  biomedieal  diagnostie  applieation,  ineluding  non-invasive  examinations  of  stranded  animals. 


RELATED  PROJECTS 

None 
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